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On the complex X-ray structure tracing the motion of Geminga 
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Abstract. A deep (100 ks) XMM-Newton observation of Geminga has shown two faint tails of diffuse X-ray emission, extend- 
ing for ~ 2' behind the pulsar, well aligned with the proper motion (PM) direction. We report here on a recent ~ 20 ks Chandra 
observation, which unveils a new structure, ~ 25" long and ~ 5" thick, starting at the pulsar position and perfectly aligned with 
the PM direction, with a surface brightness ~ 40 times higher than that of the XMM Tails. The Chandra comet-like feature has 
a remarkably hard spectrum (photon index ~ 0.9 - 1.4) and a luminosity of ~ 5.5 x 10 28 erg s -1 , comparable to the energetics 
of the larger XMM one. Geminga is thus the first neutron star to show a clear X-ray evidence of a large-scale, outer bow-shock 
as well as a short, inner cometary trail. 
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1. Introduction 

Isolated Neutron Stars (INSs) inherit high space velocities from 
their supernova explosions. Moreover, INSs are known to be ef- 
ficient particle accelerators. They power a particle wind which 
is supposed to account for the bulk of their observed rotational 
energy loss (E rot ), When the particle wind from a fast moving 
INS interacts with the surrounding Interstellar Medium (ISM), 
it gives rise to complex structures, globally named "Pulsar 
Wind Nebulae" (PWNe) where ~ 10~ 5 - 10~ 3 of the INS E rn , 
is converted into elect romagnetic radiation (see lGaensler et all 
I20041 iGaensleri I2005I for recent reviews). The study of PWNe 
may therefore give insights into the geometry and energetics 
of the particle wind and, ultimately, the configuration of the 
INS magnetosphere and the mechanisms of particle accelera- 
tion. Moreover, PWNe may probe the surrounding ISM, allow- 
ing to measure its density and ionization state. 

A basic classification of PWNe rests on the nature of the 
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external pressure confining the neutron star wind (Gaensler, 
I2005I) . For young INSs (< few 10 4 y) the pressure of the sur- 
rounding supernova ejecta plays a major role and a Crab-like 
PWN is formed. For older systems (> 10 5 y) the INS moves 
through the unperturbed ISM and the wind is confined by ram 
pressure to form a Bow-shock PWN. 

Crab-like PWNe (see ISlaneL I2005L for a review) show 
usually complex morphologies, such as tori and/or jets, typi- 
cally seen in X-rays. A remarkable axial symmetry, observed 
in several cases, is thought to trace the INS rotational axis. 
The alignement between the X-ray jets and the INS proper 
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motion directions, observed for the Crab and V ela pulsars 
llCaraveo & Mignanil 1 19991 ICaraveo et all 2001) implies an 
alignement between the rotational axis and the proper mo- 
tion of the two neutron stars, with possible important impli- 
cations Jrajheunderetanding of supernova explosion mecha- 
nisms ( lLaietalll200ll) . 



Bow-shocks (see IChatteriee & Cordesl 120021 iGaensleil 

2005, for reviews) are observed around older, less ener- 
getic INS and have a simpler, "velocity-driven" morphology. 
They are seen in H a as arc-shaped structures tracing the 
forward shock, where the neutral ISM is suddenly excited. 
Alternatively, X-ray emission (and/or radio emission on larger 
scales) is seen to trail the INS, forming a comet-like tail, due to 
synchrotron radiation from the shocked INS particles diffusing 
downstream (only in the case of PSR B 1957 +20 b oth the H„ 
and t he X-ray structures have been observed. IStappers et all 



2003). In X-rays, Bow shocks are typically fainter than Crab- 
like PWNe, favouring the detection of nearby objects. 

Proximity plays a key role in the case of Ge minga. A deep 
100 ks observation with XMM-Newton allowed lcaraveo et all 
(2003) to detect two very faint patterns, or "Tails", of diffuse 
X-ray emission trailing the neutron star, well aligned along the 
proper motion direction. 

Recently, the Geminga field was imaged by Chandra. In 
the following sections we will report on our analysis of the X- 
ray data, including a summary of the XMM-Newton results. 
We will then discuss possible physical interpretations of the 
complex morphology of the Geminga PWN. 
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Fig. 1. The left panel shows the EPIC/MOS image of the inner Geminga field. Data from MOS1 and MOS2 cameras have been 
combined. The presence of two Tails of diffuse emission, aligned with the pulsar proper motion direction (marked by the arrow), 
is apparent. The right panel shows the unfolded spectrum of the diffuse Tails (lower plot) compared t o the Geminga pulsar 
unfolded spectrum (upper plot). To describe the pulsar spectrum we used a double-component model dCaraveo et all 12004). 
encompassing a blackbody (T~ 5 x 10 5 K, green curve), representing thermal emission from the star surface, as well as a power 
law (photon index 1.7 + 0.1, blue curve), representing non-thermal emission originating in the pulsar magnetosphere. The tails' 
spectrum is well described by a hard power law (see text). 



2. The X-ray data and their analysis 

2.1. XMM-Newton observation 

Geminga has been imaged on April 4 th , 2002 for ~ 100 ks using 
the EPIC instrument. The analysis of the data collected with the 
MOS cameras yielded the detection of a symmetric structure 
consisting of two tails ~ 2 arcmin long trailing the pulsar, well 
aligned with its proper motion vector (see Fig. [2 left panel). 
The spectrum of both tails, extracted from the ~ 2 square ar- 
cmin region where the diffuse emission is resolved from the 
pulsar Point Spread Function (PSF) wings, is well described by 
a power law with an absorption consistent with that observed 
for th e pulsar itself, A^h = 1.1 X 10 20 crrT 2 (see lCaraveo et all 
120031 and references therein). Using such an Nn value, the pho- 
ton index of the 450-counts tails spectrum is 1.6 ± 0.2, simi- 
lar to the non-thermal component detected in the pulsar emis- 
sion, also shown in Fig. [0 (right panel). The observed aver- 
age surface brightness in the 0.3-5 keV range is ~ 10~ 14 erg 
cm' 2 s _1 arcmin -2 ; the total unabsorbed flux from the region 
is ~ 2.1 x 10~ 14 erg crrT 2 s _1 . At the 160 pc parallax distance, 
this translates into a 0.3-5 keV luminosity of ~ 6.5 x 10 28 erg 
s _1 , or ~ 2% of the total X-ray source luminosity in the same 
range, corresponding to a few 10~ 6 of Geminga's rotational en- 
ergy loss. 

2.2. Chandra observation 

The Chandra observation of the field of Geminga started on 
2004, February 7 at 13:02 UT and lasted for 19.9 ks. The 
Geminga pulsar was imaged on the bac k-illuminated S3 chip 
of the ACIS detector dBurke et allll997l) . 

In order to reduce the pile-up at the Geminga pulsar posi- 
tion, the 1/8 subarray mode was used. The choice of such an 
operating mode yields a time resolution of 0.7 s, but limits the 



Field Of View (FOV) to ~ 1 arcmin along the CCD readout 
direction. 

Data were retrieved through the Chandra X-ray Centre 
(CXC) Archive and were processed with the CIAO software 
v.3.2.1, using CALDB v.3.0.1, to produce calibrated "level 2" 
event lists. No periods of high background were identified, for 
a total good exposure time of 18.7 ks. In our analysis we used 
only events from the ACIS-S3 chip in the 0.45-7 keV range. 

The resulting ACIS image, zoomed on the Geminga posi- 
tion, is shown in Fig.|2](left panel) at its full resolution (pixel 
size of ~ 0.5 arcsec). In the middle panel of the same figure, 
the ACIS image, rebinned to a pixel size of ~ 2", has been 
superimposed to the XMM-Newton MOS 0.3-8 keV image. 

The Geminga pulsar is clearly detected by Chandra with 
a count rate of 7.1 + 0.2 x 10~ 2 counts s~ ! (0.45-7 keV, ex- 
traction region of 3" diameter) and a spectrum consistent 
(within t he limited sta t istics) with the results from XMM- 
Newton dCaraveo et all 120041) . The source coordinates (a = 
06 h 33 m 54.21 s , 6 = 17°46'14.2") are within ~ 0.7" from the 
values expected on the ba sis of the source absolute op tical po- 
sition and proper motion dCaraveo et allll99dll99a> . A faint 
pattern of diffuse emission, with a length of ~ 25" and a width 
of ~ 5", is also seen to trail the pulsar, perfectly aligned with 
the neutron star proper motion direction. We shall refer to such 
feature as the "Chandra Trail", while the diffuse emission de- 
tected by XMM-Newton will remain the "XMM Tails". 

The brightness profile of the Chandra Trail along the direc- 
tion of the pulsar proper motion is shown in Fig. [5] where each 
bin corresponds to the counts in a rectangular area 3" long in 
the PM direction and 7" wide perpendicular to it (see inset of 
Fig-E); the background level was evaluated from a 60" x 90" 
region North of the source region, entirely read by the same 
CCD node. In order to quantify the contribution from the bright 
pulsar in the bins immediately surrounding its centroid, an ad- 
hoc PSF was simulated using the ChaRT software, following 
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Chandra/ACIS 
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Fig. 2. The ACIS image (0.45-7 keV, 0.492" pixel size), zoomed on the Geminga pulsar position, is shown in the left panel. In 
the middle panel, the same image, rebinned to a pixel size of 2", has been superimposed on the XMM-Newton/MOS smoothed 
image. Surface brightness contours for the XMM image have been also plotted. The pulsar proper motion direction is marked by 
an arrow. The ACIS field of view is marked by the white rectangular box. The black box identifies an identical 1' X 1' region in 
all panels. In the right panel, the ACIS images in the 0.45-0.8 keV and 0.8-7 keV ranges are shown. While the Chandra Trail is 
virtually invisible in the softer band, it appears very clearly at higher energy, pointing towards a hard spectral shape. 
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Fig. 3. Brightness profile of the diffuse Chandra Trail, com- 
puted along the pulsar proper motion direction. Each bin con- 
tains the counts detected in a 7" x 3" rectangular box, as shown 
in the inset. The region represented in the inset is the same 
marked by the black box in Fig. [2 



the CXC threads and using the known spectrum of Geminga. 
The estimated profile of the point source within the histogram 
bins has been plotted in Fig.[3]as a dotted line. Although a mild 
pile-up is possibly distorting a little the actual PSF, Fig.[3]sug- 
gests that no significant diffuse emission is detected ahead of 
the pulsar. Behind the pulsar, the Chandra Trail extends along 
the proper motion direction up to a distance of ~ 25", close 
to the limit of the instrument FOV. Considering a trail region 



22" long and 7" thick extending behind Geminga, and neglect- 
ing the region within 4.5" of the pulsar, a total (background- 
subtracted) of 47 ± 7 counts (in 0.45-7 keV energy range) are 
detected. 

The bulk of the Chandra Trail counts is found above 0.8 
keV. Fig.[2](right panel) shows the very hard nature of the Trail 
emission. Although the small number of photons hampers a 
precise charachterization of the Chandra Trail emission param- 
eters, we used the maximum likeli hood metho d implemented 
in Sherpa with the CSTAT statistic dCashlll979l) to fit the spec- 
tra with an absorbed power law model, fixing the Nh to the 
value derived from X-ray fits of the Geminga pulsar spectrum 
(Nh = 1.1 x 10 2l) cm 4 ). The model describes well the data 
(the Q value - namely, the probability to observe the measured 
CSTAT statistics or higher if the model is true - is 0.65) and 
yields F = l.l*o'| (lcr). Such a range is compatible with (if 
somewhat harder) the photon index seen for the XMM Tails. 
The observed Chandra Trail flux in the 0.45-7 keV range is 
~ 2.3 x 10~ 14 erg crrT 2 s _1 . We note that the Chandra Trail 
contributes ~ 7% of the Gemin ga non-thermal flux as mea- 
sured by XMM-Newton dCaraveo et all 120041 iDe Luca et all 
2005), implying that the actual Geminga pulsed fraction above 
2 keV is slightly higher than previously reported. However, the 
correction is smaller than the statistic al error on the h igh en- 
ergy pulsed fraction quoted by IDe Luca et all J2005I) . In the 
0.3-5 keV range, the Chandra Trail has a surface brightness of 
~ 4 x 10~ 13 erg cirT 2 s _1 arcmkr 2 , i.e. ~ 40 times higher than 
the XMM Tails' one. The unabsorbed flux from the Chandra 
Trail corresponds to a luminosity of ~ 5.2 x 10 28 erg s _1 , which 
is very similar to the estimated XMM Tails' energetic. Thus, 
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the fraction of the pulsar E ro , converted into X-ray radiation 
in the arcsec Chandra Trail and in the arcmin XMM Tails is 
comparable. 

The limited FOV of the Chandra observation hampers 
a direct comparison between Chandra and XMM datasets. 
However it is worth noticing that in the immediate surrounding 
of the INS the sensitivity reached by the short Chandra obser- 
vation would not allow the detection of a diffuse emission as 
faint as the XMM Tails. Indeed, we estimate the Chandra up- 
per limit for diffuse emission at a level of ~ 10" 13 erg crrr 2 s _1 
arcmin -2 (0.45-7 keV, using a 22" x 7" region), correspond- 
ing to a flux 10 times higher than that reported for the XMM 
Tails. Thus, additional diffuse emission with a surface bright- 
ness up to 10 times the value detected by XMM-Newton could 
easily hide in the Chandra data, making it difficult to obtain a 
complete picture of the structures around Geminga. 

3. Discussion 

While the Chandra observation is both too short and with a 
FOV too limited to shed light on the nature of the XMM Tails, 
it has added yet another piece to the puzzle of Geminga. A 
bright, 25"-long trail of hard X-ray emission is now clearly 
seen behind the pulsar, well aligned with the proper motion 
direction. 

Thus, Geminga joins a small handful of pulsars sporting 
elongated, "cometary" X-ray diffuse structures, su ch as, e.g., 
the Mouse (PSR J 1747-295 8 foaensler et all 12004, the D uck 
(PSR B 1757-24. iKasni et all l200lt Idvaramadzel l2004 . the 
Black Wido w (PSR B 1957+ 20 jstappers et alll2003l) and PSR 
B 195 1+32 JLi et all l2005h which have all been interpreted 
within a pulsar bow-shock scenario. 

Although similar to the case of PSR B 195 1+32 (em- 
bedded i n the complex morphology of the SNR CTB 80, 
iLi etall [2005), the Geminga X-ray emission has a unique 
symmetry around the pulsar direction of motion. Indeed, 
the overall morphology of the diffuse features surrounding 
Geminga, with the small Chandra Trail inside the bound- 
aries of the larger, arc-shaped XMM Tails, is reminiscent 
of the composite HJX-iay nebula associated to the Black 
Widow. Such nebula is quoted as the most spectacular (and 
so far unique) confirmation of recent bow-shock theory (e.s 



Bucciantini, 2002; van der Swaluw et al., 2003; Gaensler et al., 
2004lBucciantini et allbOOsTT reflecting the double shock na- 



ture of the interaction of a fast moving neutron star with the 
ISM. An outer arc of H a radiation is seen, tracing the forward 
shock, the observed emission coming from collisionally ex- 
cited, shocked ISM; an inner, cometary X-ray feature is also 
seen, tracing the region closer to the pulsar, where the shocked 
pulsar wind emits synchrotron radiation. The case of Geminga 
show an obvious and striking difference: both the outer, arc- 
shaped, and the inner, cometary, features are seen in X-rays, 
with hard, synchrotron-like spectral shapes. 

The XMM-Newton Tails were interpreted as a Bow- 
shock, traced by very energetic (~ 10 14 eV) electrons accel- 
erated by the INS, girating and emitting synchrotron radia- 
tion in the sho c k co mpressed ambient medium magnetic field 
all 12903). Assuming a spherical pulsar wind in 



a homogeneous IS M and applin g a simple three-dimensional 
bow-shock model (Wilkin, 1996) to the XMM data it was pos- 
sible to provide an estimate of the unresolved forward stand-off 
angle 6/ w , given by the balance between the ram pressure and 
the wind pressure (0f W ~ 20"-30") and to constrain the inclina- 
tion angle of Geminga proper motion with respect to the plane 
of sky (i<30°) as well as the total ISM density (n~0.1/cm 3 ). 
Indeed, if 6/ w ~ 20" - 30" the Chandra Trail would be lo- 
cated in the unshocked pulsar wind re gion, where no stron g 
synchrotron emission is expected (e.g. Gaen sler et all 12004). 
However, owing to the large uncertainty involved in fitting the 
bow-shock model profile to the XMM diffuse feature (which is 
unresolved from the bright pulsar in its whole "forward" por- 
tion), the actual value of 6/ w could be lower and the Chandra 
Trail could still fit in the proposed scenario of pulsar wind con- 
finement. A forward stand-off distance 6/ w of ~ 10" could 
be marginally consistent with both XMM and Chandra data. 
While implying a revision of the system parameters (e.g. the 
density of the ISM would be ~ 0.5 cm -3 ), such a scenario 
would allow to interpret the Chandra Trail as the surface of 
the termination shock (i.e. the inner interface with the un- 
shocked pulsar wind cavity), where the pulsar wind particles 
are shocked and accelerated, and start emitting synchrotron ra- 
diation. The same interpretation was proposed for other ex- 
amples of cometary X-ray nebulae explained wit hin a bow- 
shock scenario, e.g. the Tongue of the Mouse JGaensler et al. , 
2004) and the nebula trailing PSR B 1757-24 jGvaramadzel 
2004), which have indeed a similar morphology. The hard spec- 
trum of the Chandra trail (F ~ 0.9 - 1 .4) is consistent, within 
the large uncertainty, with synchrotron radiation by freshly in- 
jected electrons in the slow cooling regime expected by stan- 
dard shock theory JChevaliell l2000l) . The overall shape of the 
termination shock is elongated since the ram pressure, con- 
fining the pulsar wind in the forward direction, is larger than 
the pressure in t he direction opposite to the proper motion. 
As discussed by iGaensler et alJ d2004 . this is also expected 
on the basis of detailed simulations of bqw-shqck systems (e. 



on the basis or detailed simulations or bow-shock systems (e.g. 
B uc ciantinill2002Hvan der Swaluw et all 120031 IGaensler et all 



20041 iBucciantini et allEooll) . 

Thus, within such an interpretation the Chandra Trail 
follows the surface of the pulsar wind termination shock, 
while the XMM-Newton Tails trace the larger scale, limb 
brightened, interaction of the pulsar wind with the ISM. 
Such an interaction, resulting in the formation of a bow- 
shock, could take place either through a direct contact lead- 
ing to a mixi ng between t he INS wind and the compressed 
ISM JCaraveo et all l2003l) . or through a contact discontinu- 
ity between the two shocked media, as expected by recent 
simulations (Bucciantini, 2002; van der Swaluw et al., 2003; 
IGaensler et alll2004lBucciantini et alll2005h . In the latter sce- 
nario, the XMM tails could be associated either to the shocked 
wind region (where electrons of the pulsar wind emit syn- 
chrotron radiation), or to the outer region of the shocked ISM, 
downstream of the bow-shock interface with the unperturbed 
ISM (i.e., the region which is usually associated to H a emis- 
sion in other cases of pulsar bow-shocks). Hard X-ray emission 
could be produced there via the synchrotron process by ener- 
getic electrons of the shocked pulsar wind leaking into such 
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a region, permeated by the interstellar compressed magnetic 
field, possibly higher than the magnetic field in the shocked 
wind region. 

Totally different scenarios, where the XMM Tails and the 
Chandra Trail are unrelated, could also be explored. We note 
that the phenomenology of the Chandra trail is reminiscent 
of the jet-like collimated structures seen in the case of Vela 
(Helfand et al., 2001; Pavlov et al., 2003), another example of 
nearby INS shining in high energy y-rays. In particular, the 
Geminga's Chandra Trai l can be comp ared to the "inner coun- 
terjet" of the Vela PSR dPavlov et all 120031) in terms of their 
absolute projected dimensions (~ 5 x 10 16 cm), spectral shape 
(photon index ~1.2) and efficiency (Lx~l0^ b E mt ). Although 
the complex features of the Vela PWN still lack a consistent ex- 
planation, and their 3-D morpho logy and orientation is far from 
understood dPavlovet alll2003l) . they clearly show that the as- 
sumption of an isotropic wind outflow is not adequate. As in 
the case of Vela, Geminga's diffuse features are characterized 
by a definite symmetry along the proper motion direction. 

Before drawing conclusions on the morphology of the 
structures trailing Geminga, we would like to underline that 
the combined Chandra/XMM-Newton image shown in Fig. [2] 
however tantalizing, does not provide a complete and unbiased 
picture of the INS surroundings. We see bright emission next 
to the INS and faint emission far from it but, owing to instru- 
ment characteristics and/or observing time limitations, Fig. [5] 
lacks sensitivity in the intermediate flux region where addi- 
tional structures may be present. Only a deep Chandra observa- 
tion could yield a complete picture of the system morphology, 
especially in its forward portion, the piece of evidence needed 
to ultimately understand its physics. 
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